Abstract 


Presentation 


Table  of  Contents 


Next  Abstract 


LINEAR  FEATURE  DETECTION  IN  S AR  IMAGES 

Ph.  COURMONTAGNE 
L2MP  UMR  CNRS  6137  /  ISEM 


CENTRE  NATIONAL 
DE  LA  RECHERCHE 
SCIENTIFIQUE 


UMR  6137 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
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LINEAR  FEATURE  DETECTION  IN  SAR  IMAGES 


INTRODUCTION 


O  Ship  wakes  in  Synthetic  Aperture  Radar  (SAR)  images 
O  Detection  using  Radon  transform  (DEANS,  1983) 


O  SAR  images  affected  by  speckle 

^  wedding  between  the  Radon  transform  and  a  filtering  method 
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LINEAR  FEATURE  DETECTION  IN  SAR  IMAGES 


CONTENTS 

O  Radon  transform 

^  discrete  form 

O  Stochastic  matched  filtering  method 

O  Interpolation-filtering  method 

^  theory  and  subimage  processing 

O  Experimental  results  on  SAR  images 

^  comparison  with  the  classical  approach 
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THE  RADON  TRANSFORM 

O  Consider  an  image  /,  with  dimensions  M  x  M 
The  Radon  transform  I  of  this  image  is: 

M/2 

I{x0,9)=  I{x0cos6  —  y0sinO,X0sm9 t/ffCOsO) 

y0=-M/2 

where  {xg  ,  yo)  e  Z  and  0  C  [0;  tt] 

O  Problems  :  -  new  pixel  coordinates  in  3?  ^ 


^  Computation  of  the  new  pixel  values 


-  pixels  localized  at  the  edge  in  ^  I  (surrounding  points) 

^  Edge  effect  size  of  the  image  in  ^  x  ^ 

\/2 
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STOCHASTIC  MATCHED  FILTERING  METHOD 


OConsider  the  stationary  noise -corrupted  signal  defined  overD : 

Z{x,y)  =  S{x,y)  +  B{x,y) 

where  signal  S{x,  y)  and  noise  B{x,  y)  are  assumed  to  be  independent. 
O Observed  signal  expansion: 

Q 

Z{x,y)  =  '^z„'^n{x,y) 


n=\ 


where:  / 


:  deterministic  and  linearly  independent  basis  functions 
Q  :  number  of  basis  functions  retained  for  the  expansion 


uncorrelated  random  variables  defined  by  2:^  =  /  /  Z {x,  y)dxdy 

<  J  J  D 
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LINEAR  FEATURE  DETECTION  IN  SAR  IMAGES 


O Determination  of  functions  v) 

"^Making  sure  Zn  are  uncorrelated 

'^Optimization  of  the  signal  to  noise  ratio  K  expressed  as  a  Rayleigh  quotient 
=>  K  will  be  maximized  if  v)  are  solutions  of: 


JJ  rss{x-x\tj-ij')^n{x\y')dxdij=Xn  jj  TBB{x-x,y-tJ)^n{x\y)dx'dy 


where  F and  Tbb  are  the  covariances  of  the  signal  and  of  the  noise 


O Determination  of  functions  y) 


Ij)  =  JJ  rBB(x  -  x',  y  -  tj)^n{x\  ij)dx'dy' 


O Signal  to  noise  ratio  of  the  component  of  Z(x,  ^) :  ^  An 
=>  number  Q  of  basis  functions  such  as  Ag  >  1 
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INTERPOLATION-FILTERING  METHOD 


'^Observed  signal  expansion  and  restoration  using  interpolated  funetions  v) 


O  Defaults:  heavy  CPU  budget  and  memory  problems 

=>  new  formulation  using  the  diserete  eosine  transform 

- L2MP  Laboratory  — 
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O  Z{x,  y)  interpolation-filtering  using  DCT,  with  D  =  [— T;  Tf’ 


^n{x,  y)  <^= 

Nf  Nf  Nf  Nf 

X  X  X  X 

k=0  1=0  k=0  1=0 

'^n{x,y)  <^= 

Nf  Nf 

A:=0  1=0 

Nf  Nf 

^  X  X 

p=0  q=0 

Z{x,y)  <^= 

Q 

=>  ^  ZnPk,l 

n=l 

Functions 

DCT  coefficients 

^n{x,y) 

Z{x,  y) 

d 

^p,q 

Z(x,y) 

dk,i 

Ui- 

Thb 

O  Signal  of  interest  restoration  in 

Nf  Nf 

Ze{xe,ije)  =  XI XI 

k=Q  /=0 


Txk[xe  —  T) 
2T 


cos 


TTljye  -  T) 
2T 


9th  ASAP  Workshop  13/14  March  2001  -  7/17 


L2MP  Laboratory  — 


LINEAR  FEATURE  DETECTION  IN  SAR  IMAGES 


'll? 

kU 


Subimage  processing 

OAssumption:  Z[x^y)  is  stationary 
=>  necessity  of  a  subimage  processing 


OEdge  effect 


overlapping  between  adjacent  subimages: 


N  ,  ,  N  sin  6 

d  = -  and  h  =  - - 

cos  6>  +  sin  (9  sin  6  cos  0 


O  Smoothing  effect  when  Q  is  the  same  for  the  whole  image 

=>  minimization  of  the  mean  square  error  to  find  Q  for  each  subimage 

.  Q  Nf  Nf 


^  +  7  X!  (c'l  - 


iPh) 


k=0  1=0 
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EXPERIMENTAL  RESULTS 

OSIR-C/X-SAR  image 
'^moving  ship 

'^dark  turbulent  wake  (17  miles) 


0256  gray  levels  (0:  black,  255:  white) 


Olmage  size:  698  X  698 


O  Variation  coefficient:  0.277 


-300  -200  -100 


100  200  300 


ODifficulty:  dark  patches  in  the  upper  right  corresponding  to  smooth  area  of  low  wind 
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Signal  and  noise  auto -correlation  functions 

Oa  priori  knowledge  of  the  signal  and  the  noise  auto-correlation  functions 
=>  determination  of  normalized  auto-correlation  models 


t  .2  -  1  2  - . 


Signal  auto-correlation  function  Noise  auto-correlation  function 
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Interpolation-filtering  of  the  SIR-C/X-SAR  image 


O Rotation  angle:  35*^ 

ONumber  Q  of  basis  functions: 

^  approximately  13  for  the  wake 
^  near  1  for  the  rest  of  the  image 

Olmage  size:  973  X  973 

O  Variation  coefficient:  0.016 


SIR-C/X-SAR  image  in  reference  system  §^35 
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Radon  domain  displayed  as  image 


-400  -300  -200  -100  0  'OO  200  300  40O 


Projections 
Our  processing 


-<oo  -300  -200  -‘00  0  '  00  2X  an  «xi 


Projections 
Classical  processing 


2 


^  Trough  corresponding  to  the  wake  localized  near  51^ 
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Amplitude 


LINEAR  FEATURE  DETECTION  IN  SAR  IMAGES 


Radon  domain  displayed  as  graph 


ambiguity  about  the  presence  of  a  ship  wake  with  our  processing 

=>  detection  by  simply  using  a  threshold  contrarily  to  classical  approach 
- L2MP  Laboratory 
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Reconstructed  image  using  the  transform  domain 


'^Inverse  Radon  transform  applied  to  the  transform  domain  raised  to  the  power  of  3 


Our  processing 
Variation  coefficient:  0.002 


Classical  processing 
Variation  coefficient:  0.145 
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ANOTHER  SAR  IMAGE 


-150  -100  -50  0  50  100  150 


OERS  SAR  image 

'^ship  pixels  replaced  by  mean  background 
"^dark  turbulent  wake 


0256  gray  levels  (0:  black,  255:  white) 


Olmage  size:  325  x  325 


ODifficulty:  dark  patch  near  the  wake  (slick  oil) 
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Radon  domain  displayed  as  image 


Projections 
Our  processing 


-200  -150  -'00  -50  0  50  100  -50  200 

Projections 


Classical  processing 


^  Trough  corresponding  to  the  wake  localized  near  11 5^ 
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CONCLUSION 

O  New  processing  for  ship  wakes  detection 
^  SAR  image  Radon  transform 

^  original  contribution:  taking  into  account  the  speckle  for  image  interpolation 

O  Better  detection  with  lower  probability  of  false  alarm  or  no  detection 

O  Future  work:  extension  of  this  method  to  the  localized  Radon  transform 
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